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ABSTRACT 

The purpose of this letter is threefold : (i) to derive, in the framework of a new parametrization, some 
compact formulas of energy averages for the electrostatic interaction within an n£ N configuration, (ii) to describe 
a new generating function for obtaining the number of states with a given spin angular momentum in an nt N 
configuration, and (iii) to report some apparently new sum rules (actually a by-product of (i)) for SU(2) D U(l) 
coupling coefficients. 
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1. Introduction 

In the theory of complex spectra [1-3], one- and two-body Hamiltonians invariant under the group 50(3) 
and symmetric in the spin and orbital parts can be written in the form [4-5] 

V = J2Y, D[(kik 2 )k s (k 3 k4)k L } {{ M ( fcl )W® u ( fc ^(j)} (fcs) «'{u ( ' C3) (i)«'u ( ' C4) (j)} ( ' Ct) }o 0) - (1) 

i^j all k 

In eq. (1), stands for a Racah unit tensor of rank k. Further, the 50(3)-invariant operator { }q ^ results 
from the coupling of the tensor product {u^- kl \i) <g>u( k2 \j)}( ks ' > acting on the spin part with the tensor product 
{u^ k3 \i) (g>u(' C4 - ) (j)}( fcl '- ) acting on the orbital part. Finally, the D[ ] parameters in eq. (1) are radial parameters 
(depending on the radial wavefunctions involved, e.g., R n e(r), R n i^(r), etc.) which are generally taken as 
phenomenological parameters. 

The operator (1) can be considered as a particular case of the G- invariant Hamiltonian introduced in 
[5] for describing optical and magnetic properties of a partly-filled shell ion in a crystalline environment with 
symmetry G. Equation (1) corresponds to G = 50(3) : to obtain (1) from ref. [5], it is enough to put k = 
and to replace a or a r 7o by q = 0. 

We shall be concerned in this work with the (spin-independent) Coulomb interaction which is obtained 
from (1) by taking k\ = k 2 = and fc 3 = k± = k. In this case, the parameters J5[(00)0(fc/c)0] are proportional to 
the Slater parameters F^ k \ We shall restrict ourselves to the action of the Coulomb interaction V{k\ = k 2 = 
0) within an nl N configuration but, for the purpose of forthcoming generalizations, we shall consider Slater 
parameters of the type (£,£') — {£,£'; £,£'). The F^ {£,£') parametrization corresponds to a multipolar 
expansion of the electrostatic interaction V(k\ — k 2 — 0) and turns out to be especially adapted to the chain 
50(3) D 50(2). 

There are several other parametrizations besides the F^ k \£,l') parametrization. (For instance, the E k 
parametrization [3] is well-known for nf N configurations.) An alternative parametrization, referred to here as 
the the £ x (£, £') parametrization, was introduced in ref. [6]. This parametrization was obtained by transposing to 
atomic and nuclear spectroscopy a parametrization, namely, the Angular Overlap Model parametrization [7,8], 
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used in the spectroscopy of partly-filled shell ions in crystalline environments. The £ x (£,£') parametrization 
can be defined by the passage formula 

or the reverse formula 

£\ZJ>) = {-lf^W+Wt r +V E(o o)(-A I) F(fc) ^ £ ')' ( 3 ) 
which generalizes eq. (2) of ref. [6]. 

It is one of the aims of this letter to show (in section 3) that various energy averages for nl N configurations 
assume a particularly simple form when expressed in the £ x {£,£') parametrization. To obtain energy averages, 
one needs to know the number of states, having well-defined qualifications, in the n£ N configuration and a new 
way of denumbering states, originally introduced in ref. [9], is further developed in section 2. Two new sum rules 
for 3 — jm symbols, which are at the root of the derivation of two energy averages in section 3, are relegated to 
an appendix. Finally, some concluding remarks are given in section 4. 

2. Denumbering states 

Let us consider a system of N fermions in a shell n£. As is well-known, the resulting configuration n£ N 
has ( 4 ^" 2 ) totally anti-symmetric state vectors. Among these ( 4 ^ 2 ) states, let Hi(N, S) be the number of states 
having a given total spin S. We devote this section to the calculation of He(N, S). 

From ref. [9] , we know that the function 

1/2 I 

F(x,y,z)= [] [] (l + ^ m ^ mf ) (4) 

m s = -l/2 m e = -£ 

is the generating function for the number Fe(N, Ms, Ml) of states in n£ N with z-components of the total spin 
and orbital angular momenta equal to Ms and Ml, respectively. The number Fg(N, Ms, Ml) is obtained by 
expanding F(x,y,z) as 

F(x,y,z)= Yl F i (N,M s ,M L )z N y Ms x M K (5) 

N , Ms, M L 
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By setting x = 1, we obtain a generating function for 

G e (N, M S ) = ]T F e (N, M s , M L ) (6) 

M L 

which is the number of states with a definite value of M$. The latter generating function can be written in two 
equivalent forms, namely : 

r i 21+1 

(l + zy- 1 ' 2 ) 21 * 1 (1 + zy 1 ' 2 )"* 1 or [l + z (y^ 2 + y 1 ' 2 ) + z 2 ] . (7) 
This leads to the following expression 

WM S )=(| + ^)(| + (8) 

or alternatively 

GtW M s) = (M+ l)l g a(2l+1 _ N + m _ i _ M3m _ i + Us)r (•) 

In particular, the compact form (8) is very simple to handle. Finally, the number Ht(N,S) of states of the 
configuration nl N with total spin S is simply obtained by combining 

H e (N,S) = Ge(N,S)-G e (N,S + l) for 5 < y, 

AT AT AT 

H t (N,-)=G t (N,-) for 5= y (10) 

with cq. (8). 

As an illustration, we consider the configuration nf 6 . From eq. (8), the numbers G 3 (6, Ms) are found to 
be 7, 147, 735 and 1225 for M s = 3, 2, 1 and 0, respectively. Furthermore from eq. (10), the numbers H 3 (6,S) 
are easily seen to be 7, 140, 588 and 490 for 5 = 3,2,1 and 0, respectively. As a check, we verify that the total 
number of states is 3003 = (g 4 ) . 

3. Averages for nl N configurations 

3.1. Average interaction energy 



The average interaction energy E av (£, £) for an arbitrary nl N configuration in spherical symmetry was 
derived by Shortley [2] and further discussed by Slater in his book [1] (see also the nice book by Condon and 
Odaba§i [2]). The expression for E av (£,£) is known in terms of the Slater parameters F k {£,£). The formula 
for E av (£,£) in the F k (£,£) parametrization [1,2] does not exhibit any remarkable peculiarity. It is a simple 
matter of calculation (by means of Wigner-Racah calculus) to convert the expression for E av (£, £) in the £ x (£, £) 
parametrization. This yields 

It is to be noted that, from a fitting procedure viewpoint, cq. (11) involves two parameters (£ a = £° and 
^2x =1 £ x ). This inclines us to introduce the linear combinations 

s = \jz£\ v = ih {£rr - s) - (12) 

A=l 

Thus, cq. (11) can be rewritten as 

Em(lJ)= ^i^ s+ ^± v y (13) 

in terms of the non-independent parameters S and V. 

In the special case where all the parameters £ A are taken to be equal, say to a test value £ (i.e., S = £ 
and V = 0), equations (11) and (13) lead to 

E av (£,£)= N{N ~ ^ £. (14) 

Such a kind of result is especially important for the purpose of checking electrostatic interaction matrices. 
Indeed, in the case where £ x = £ for any A, it can be shown that all energy levels have the same value, viz., 
[N(N— l)/2]£ ; then, the spectrum of V(ki = k 2 = 0) is maximally degenerate and (14) is a simple consequence 
of this maximal degeneracy. 

3.2. Other average energies 



N N 
_(_ + !) 



We now turn our attention to the average energy 2S+1 E(£, £) over all the states of the configuration n£ N 
corresponding to a given total spin S. From extensive calculations of the electrostatic energy levels for the 
electronic configurations np N , nd N and nf N , we empirically discovered that E av (£,£) is given by 

V. (15) 

(All matrix elements of the Coulomb interaction V(k\ = k 2 = 0) for the configurations np N , nd N and nf N 
listed in rcf. [10] , in the E k (Racah) parametrization for I = f and in the (£, £') (Slater) parametrization for 
£ = d and p, have been transcribed in the £ x (£, £') parametrization with the help of the algebraic and symbolic 
programming system REDUCE. The complete electrostatic energy matrices for the electronic configurations 
np N , nd N and nf N in the £ x {£,£') parametrization are presently under preparation for distribution to the 
interested readers.) The case S — -j is of special interest ; in this case, which corresponds to the highest 
multiplicity term, cq. (15) simply reduces to 

^E av (£,£) = ^fVs, (16) 

so that N+1 E av (£,£) does not depend on the parameter £" . 

Here again, we note that in the special case where £ A = £ for any A, we obtain from (15) 

■^ Em{U) = m r i> £ , (17) 

a result to be compared with eq. (14). 

The proof of eq. (15) can be achieved by constructing all determinental wavefunctions with a given value 
Ms of the z-componcnt of the total spin angular momentum and by calculating the sum of the energies of 
all these determinental wavefunctions. The complete proof shall be reported elsewhere. The proof of (15) for 
N > 2 can be also obtained, in principle, as an extension of the one for N — 2 and we now derive eq. (15) for 
an n£ 2 configuration. 

We start from the relation (3) of ref. [6] giving the energy of the term 25+1 L of the configuration n£ 2 . 
Such a relation can be rewritten as 

t 2 

^L={2£+1) ]T ( l [ L x ) £\£,£). (18) 
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Therefore, for the configuration under consideration, we have 

1 

El. + 1) 



2S +^E av {£,£) = 1 V(2^ + f) 25 + 1 ^, (19) 



where the sums over are to be performed on even values of -L, — 0(2)(2£), or on odd values of _L, 
L n = 1(2) (21 — 1), according to whether 5 is (singlet states) or 1 (triplet states). (As usual, the notation 
i = a(b)c means that i takes the values a, a + b, a + 26, . . ., a + [^j^Jfe.) By combining eqs. (18) and (19), we get 

We -M° E^Ti) £ £HU) S (2i ' + 1) (S a m) 2 ' < 2 °> 

Although there is a well-known formula for expressing the last sum in (20) when El. is replaced by El w i tn 
L = 0(1)(2£), to the best of our knowledge there is no formula in the literature for calculating the last sum in 
(20) for the two distinct cases where = Q(2)(2£) and L n = 1(2)(2£ — 1). By using the formula (30) derived 
in the appendix, the sum El^m m (20) can be calculated and we finally arrive at 

1 - t" (£,£) + 

\=i 



X E av (£,£) 



(21) 



'+1 

for the singlet states and 

1 £ 

3 Eav (e,£) = - fA (^) ( 22 ) 

for the triplet states. It is immediate to check that eqs. (21) and (22) are particular cases of the general formula 
(15) corresponding to (N = 2, 5 = 0) and (N = 2, 5 — 1), respectively. 



The consistency of (13) and (15) requires that 

E s (25+l)g,(iV,5) | [f(f + l) -5(5 + 1)] _ £+1 N(N-l) 
Es(2S + l) H e (N,S) U + l 2 



(23) 



~ 2 3N / N-l\ 
< S > = -T i-TTTT . (2 4 ) 



from which we easily deduce 

<r 5 2 ~> = - 

4 V 4£ + 1 , 

where < S 2 > is the average, over all the states of the configuration n£ N , of the square S 2 of the total spin 
angular momentum. The formula so-obtained for < S 2 > agrees with the one it is possible to derive from 

£ 2 = Z s (2S + l)H t (N,S) 5(5 + 1) 
/ U + 2 x 
N 



by using the explicit expression for Hi{N, S) given in eqs. (8) and (10). 



4. Concluding remarks 

In this paper we have concentrated on electrostatic energy averages for nt N configurations in a new 
parametrization, viz., the £ x (£,£') parametrization defined by (2) and (3). It is to be emphasized that the ob- 
tained averages (13) and (15) depend only on two parameters (<S and V). This result and the fact, already noted 
in ref. [6], that the term energies for the configuration n£ 2 assume a very simple form in the £ x (£, £') parametriza- 
tion, are two indications that a hidden symmetry is probably inherent to the £ x (£,£') parametrization. In this 
respect, it would be interesting to find a group theoretical interpretation for the £ x (£,£') parametrization. 

In addition to the well-known interest, mentioned in refs. [1] and [2], of energy averages, it is to be pointed 
out that compact formulas, like (13) and (15), for such averages constitute useful means for checking energy 
matrices. Furthermore, eq. (15) suggests a strong version of Hund's rule, according to which the energy average 
2S+1 E av (£, £) over all the states having a given total spin S decreases (linearly) in S(S + 1) upon increasing S. 
This statement depends on the fact that the coefficient of S(S + 1) in (15) should be negative. It would be of 
interest to examine the effect that an independent optimization of the radial wavefunction for the average energy 
corresponding to each spin S will have on the functional form of the dependence of this energy on S(S + 1) (see 
also ref. [11]). 

The result (15) concerns the electrostatic interaction, in the £ x (£,£') parametrization, within an nl N con- 
figuration in spherical symmetry (SL coupling). This suggests several possible extensions of some of the results 
contained in the present paper. In particular, it would be appealing to extend the £ x (£,£') parametrization 
to other interactions (e.g., to the general interaction described by (1)) and to other configurations in arbitrary 
symmetry (e.g., to atomic and nuclear configurations with several open shells in spherical symmetry or to the 
molecular configuration in cubical symmetry). Along this vein, it is to be mentioned that electro- 

static energy averages have been derived, in the (£,£') parametrization, for j N configurations in spherical 
symmetry (jj coupling) [12]. 
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To derive (15) we used (10) which furnishes a new way for obtaining the number of states in n£ N with a 
given spin S in the case of spherical symmetry (G = 50(3)). It would be worthwhile to extend (10), and more 
geneally the generating function (4), to the case of an arbitrary point symmmetry group G. It would be also 
very useful to extend (4) to other systems (e.g., quark systems) than systems involving electrons by introducing 
in (4) additional degrees of freedom (e.g., isospin, flavor and color). 

The energy averages (21) and (22) actually were obtained by inspection of tables giving term energies for 
n£ N in the £ x {£,£') parametrization. As is often the case in spectroscopy, some regurality in energy formulas is 
the signature of special relations between coupling and/or recoupling coefficients. In this regard, we showed that 
the compact formulas (21) and (22) are a direct consequence of the sum rules for Clebsch-Gordan coefficients 
derived in the appendix. 
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Appendix. Sum rules for 3 jm symbols 

The aim of this appendix is to report two apparently new sum rules for SU (2) D U(l) coupling coefficients. 
Equations (21) and (22) of the main body of this paper are simple consequences of these sum rules. 

Let us consider the decomposition 

W®W = {W®(m+0[W®W]-, 

{(*)®(m+= (ie), [V)®(l)]-= (La), (26) 

L„=0(2)(2^) L =l(2)(2£-1) 

into a symmetrized part { } + and an anti-symmetrized part [ ]_, of the direct product (£)<&(£) of two irreducible 
representation classes (£) of SU{2). Such a decomposition can be transcribed in terms of basis vectors, acting 
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on two spaces (1 and 2) of constant angular momentum £, as 

\£m) 1 ® \£m') 2 + ir |&n')i ® |*m) 2 = 2 ^ 1^ 

(^toto'I-L^M), (27) 



where 

L^ = L e = 0(2)(2f) for tt = +1 and L V = L = 1(2) (2^ - 1) for tt = -1. (28) 
By taking the scalar product of (27) with £/x)i ® |^/x')2, we get the sum rules 

5(m, //) + tt 5(m', /z) 5(m, fi?) = 2 ^ (« [ip!\L^M) {U mm'\L^M) (29) 

or, equivalently, in terms of 3 — jm symbols 

E( 2L - + 1 )( L m I i')(-M t J)4^^( m >')+^'^KA')]. (30) 
The particular case (to = /i = 0, m' = // = A) leads to (21) and (22) for 7r = +1 and 7r = —1, respectively 
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